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Abstract 

The magnetic interaction of a pair of bound magnetic polarons (BMP) in di- 
luted magnetic semiconductors (DMS) is analyzed via a generalized Hubbard- 
type Hamiltonian for two carriers in the presence of effective magnetic fields 
arising from the magnetic polarization of their respective polarons. For the 
case where the magnetic fields at the two sites have equal magnitude but are 
allowed to have arbitrary directions, it is shown that the energy of the two 
polarons is minimized for a ferromagnetic configuration of the carrier spins 
(in contrast to the case of hydrogenic centers in nonmagnetic semiconductors) 
if polaron fields are strong enough. A modified Heisenberg-type Hamiltonian 
is constructed to describe the low energy states of the resulting system. 

I. INTRODUCTION 

Shallow impurities in doped semiconductors0 can be described in terms of a hydro- 
genic Hamiltonian with an effective mass given by the band mass, and a Coulomb potential 
screened by the dielectric constant of the host semiconductor. While the "effective mass" 
equation is slightly more complicated for donors in indirect band-gap semiconductors, and for 
acceptors a matrix version of a generalized Hydrogen problem is obtained, these differences 
change details, but not the basic physics. Therefore, the hydrogenic model is a useful guide 
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for studying and understanding interactions between impurities in doped semiconductorsatl . 
At low concentrations the interactions between impurity centers can be modelled in terms of 
an exchange Hamiltonian involving, as the dominant term@, pairwise Heisenberg exchange 
corresponding to the Hydrogen molecule problem!, where the exchange interactions are 
known to be anti-ferromagnetic at all distances!. 

In diluted magnetic semiconductors (DMS), a small fraction of the nonmagnetic ions 
that form the lattice are replaced by magnetic ions such as Mn or Fe. Several features (such 
as variable band gaps, optical response, spin polarized transport, as well as the unusual 
magnetic behaviour analyzed in this paper) have turned DMS into a topic of considerable 
interest during recent years. In the low doping regime (i.e. carrier density below the Mott 
transition), the susceptibility (i.e., dM/dH) vs. magnetic field (H) curve of such a DMS 
displays a curious double-step feature!. To understand the properties of DMS with dopants, 
it is not only necessary to understand the direct interactions between the dopants, but also 
the interactions with the magnetic ions which by themselves contain low lying degrees of 
freedom. When the magnetic ions are dilute, their direct interactions are unimportant. 
Thus, for example, the problem of a single shallow impurity in a DMS is well described 
in terms of an exchange interaction between the bound carrier (electron or hole) and the 
magnetic ion, and is known from extensive studies to lead to the formation of a bound 
magnetic polaron (BMP)!. The spins of the magnetic impurity ions within one effective 
Bohr radius of a dopant interact via a sizable exchange with the carrier, thus becoming 
aligned and forming large-spin polarons. The polarons align with an external magnetic field 
before the individual magnetic ions do, thus giving rise to the two-step susceptibility curve. 
By analyzing the step due to BMPs (fitting it to a Curie- Weiss form), a ferromagnetic 
interraction between the polarons can be deduced. However this result seems puzzling: in 
non-magnetic semiconductors carrier virtual hopping invariably yields anti-ferromagnetismll 



that a ferromagnetic interraction between the polarons can be obtained if one considers the 
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overlap of two polarons formed around two dopants. The two carriers compete over the 
spins in the overlapping region. For certain values of the model parameters the indirect 
carrier-magnetic ion-carrier interaction becomes stronger than the direct carrier exchange, 
and the polarons align resulting in ferromagnetism. 

The current work approaches the problem from a different perspective. The polarons are 
considered as non-overlapping, and their effect on the dopant atom is taken into account 
through a local magnetic field h. A system of two such polarons is analyzed via a generalized 
Hubbard-type Hamiltonian, where hopping (matrix element t) and Coulomb interaction 
(energy U) are turned on. Several cases are considered (dopants with a single bound state 
and with several bound excited states). Ferromagnetic behaviour is seen to emerge when the 
carrier is allowed to hop between the ground state of one dopant atom and excited states 
of the other dopant. Numerical work supports the conclusion that such a ferromagnetic 
interaction is indeed possible in realistic conditions. We then discuss the applicability of a 
Heisenberg-type model for two interacting polarons. In the moderately high field domain 
t <C h <C U, an effective Heisenberg Hamiltonian is found which contains a mixing of the 
magnetic fields at the two sites. 



II. THE MODEL AND THE APPROACH 

Our model consists of two identical atoms, with several bound states, in arbitrary local 
magnetic fields. We allow for hopping of the carriers between the two atoms. The local 
magnetic fields represent the exchange fields due to the magnetic ions at each site. As the 
number of magnetic ions around each site is large, and the doping is considered uniform, we 
assume that the magnitudes of the magnetic fields at the two sites are equal. However, the 
directions of the two fields are allowed to be arbitrary. Thus, the Hamiltonian we study has 
the general form: 

h = -m + vi) - hHnys, - hh(r 2 ) -s 2 -^-^-^-^ + ^ (i) 
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where h is the magnitude of the field, h(r) is the (arbitrary) direction, a and b are the labels 
for the two Hydrogenic centers and 1 and 2 are the labels for the two electrons. 

If we consider a Hubbard-like approximation! with one energy level per impurity site 
and no magnetic fields, the Hamiltonian becomes (in second quantized form) : 

H = J2 ( C ( n «T + n «l) + Ull^llai) + t( C L C bs + <4 C as) (2) 

a=a,b s=T,| 

where a, b are labels for the two impurity sites, c a f is the annihilation operator for the state 
on impurity a with up-spin, n a j is the occupation number, = cLc a f etc. 

If we also introduce arbitrary number of energy levels on each impurity atom, the Hub- 
bard Hamiltonian (Q) turns into : 

H Si,a,s ^ia^ias 4~ Sj<j,si,S2,a ^ij Hj QSl llj aS2 4" 

(3) 

where a G {a, b} indexes the impurity sites; i,j the atomic levels on each impurity; Si, S2 € 
{|,|} the spin degree of freedom; e« the energy of level i; and Uij the Coulomb interaction 
energy of electrons in states i and j on the same impurity atom. 

Finally, if at each site we consider the arbitrary magnetic fields h a , hf,, and we quantize 
spin along the axes of the local magnetic fields (i.e. c\ a + creates an electron in the z-th state 
on impurity a with spin parallel to h a ), the Hubbard Hamiltonian becomes : 

H Si,a,s Cia^-ias 4~ Si<j',si,S2,a ^ij^Haai^-jass 4" Si, a U iiHi a ^\\i a y-\- 

(4) 

4" Si,a ^a( n iaT n ia|) 4" St,j,si,s 2 ^«si j's 2 ( C j6si C ias 2 "I" C ias 2 C jfrsi ) 

We caution the reader that in this case the transition matrix elements t isi j S2 become depen- 
dent on the angle between the two magnetic fields h a , h^. We will discuss the relationship 
between Eqs. |]||,@ and Eq. [I] in more detail in the concluding section. 

Several models of increasing complexity were considered : atoms with a single bound 
state and without magnetic fields (sec. [Ill A|) , atoms with a single bound state in arbitrary 



magnetic fields (sec. |1 [ [B|) , atoms with several excited states in arbitrary magnetic fields 
(sec. |1V A| ). The ground state of the two-center system is shown to undergo a transition 
from an antiferromagnetic state (singlet) to a fully ferromagnetic (triplet) configuration with 
the increase of the effective polaron magnetic field (sec. [IV B| ). The results for a regular 
Heisenberg Hamiltonian where the two spins are in arbitrary fixed fields are calculated as 
well and compared with those derived from our model (sec. [Ill C[) . Finally we find a modified 
Heisenberg type Hamiltonian that agrees with our model in the moderately high field regime 
(sec. 0). 

III. REGULAR HUBBARD MODEL 

A. Regular Hubbard model in zero field 

The Hubbard model of the hydrogen molecule! (see Eq. Q) consists of two hydrogenic 
(one-electron) centers, each with one single electron bound state of energy e. Electrons 
are allowed to hop between the two sites, subject to the restrictions imposed by the Pauli 
principle with a hopping matrix element t. Each center also has one two-electron state, with 
energy 2e + U, where U represents the interaction energy between two electrons on the same 
atom. 

In this case we can decouple the spin-1 and spin-0 subspaces (since there is no connection 
between states of different spin). It turns out that the spectrum of eigenvalues is 2e (triple 
degenerate) for spin-1 and 2e — jj-, 2e + U and 2e + U + ^ for spin-0. The ground state 
therefore has spin {i.e. the interaction between the electrons can be thought of as anti- 
ferromagnetic). 

B. Regular Hubbard model with arbitrary fields 

If two arbitrary fixed fields (of strengths h a and and making an angle O) are applied 
at the two sites, the analysis becomes more complicated (see Eq. |], and consider that there 
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is a single level, % = 1, on each impurity). We can quantize the spins along the axes of the 
fields, and we can pick as a basis: 

{ata|,aT&T; a T&l) a l&1\ a l&l)&T&l} 
Each element in the basis is anti-symmetrized, for example 

a t b 1= -^(|Va(zi) T)W*3) I) - Mm) DWaM T» 

This basis obviously is not formed of eigenstates of the total spin, but only of the z- 

component of the spin. 

The new Hamiltonian matrix (ignoring the Bohr magneton, for simplicity) is shown 

in Eq. || below : 
/ 



2e + U 



it sin — 



t cos ® —t cos y 



-it sin | 2e + 



i cos 



e 



-t cos 



e 



2e+ ^ 









it sin ~ 











it sin ® 








t cos ® 



2e - 



hg-hb 



-it sin ~ 













-it sin ^ 



t cos 
-t cos 



it sin ^ 



-t cos y 



-itsinf 2e + C/ 



(5) 



We next make the simplifying assumption that h a = = h. For the case of DMS, since 
the magnetic ion distribution is random, this can be justified if each polaron has several (N) 
magnetic ions producing the exchange field on the carrier, so \h a — h b \ = h/y/N -C h. In this 
case we obtain two pairs of degenerate states {a ] a [,b \ b [} and {a t b I, a I b t}. By 
making a 45 degree rotation within each of the degenerate subspaces, and by multiplying 
some of the basis vectors by i when necessary, we single out two of the eigenvalues (2e and 
2e + U), the rest of the matrix having the simpler form in Eq. |6| below : 



2e + h 






2e 









y/2t sin 



— 2t cos 



2e - h \/2t sin | 



\/2tsinf -2tcosf ^sinf 2e + U 
6 



(6) 



This Hamiltonian can be solved by perturbation theory (PT). In the limit h — > the 
upper three states become degenerate and the problem needs to be handled by degenerate 
perturbation theory. We will not investigate this limit any further. In the high field limit 
however, the magnetic field removes the degeneracy and we can obtain the eigenvalues to 
second order by regular PT. Thus we obtain for the lowest eigenvalues : 

fa+h-^ffi 2e 2e - 2e-h- 2 -^^) (7) 



C. Heisenberg Hamiltonian with arbitrary fields 



By solving the same problem (two atoms in fixed external fields) using a Heisenberg 
Hamiltonian, and comparing the eignevalues with the results obtained above in Eq. [7], one 
can see how the effective exchange parameter in the Heisenberg Hamiltonian is affected by 
external magnetic fields. We start with0 : 



H H i = h a ■ s a + h b ■ s b + Js a ■ s b 



(8) 



We can again quantize the spins along the axes of the two fields, and work in the basis 



{a t b |, a t b |,a j b |, a J, b j} 



which yields the Hamiltonian matrix : 
/ 



H 



HI 



—h + ^ cos 6 
-|sin9 
| sin 6 



—j sin0 



sin 9 



-£(l-cos9) ' 



-^cos9 4(1 + cos 9) -^sin9 



4 ~ 4 

^(l + cos9) -^cos9 



-~(l-cos9) — ~ sin 9 4 sin 9 



^sin9 
h + j cos 9 



(9) 



4 ^ w„ ^ j 4 ^„ w 4 

After we do a rotation by 45 degrees in the {a ] b [,a [ b ]} subspace we can apply 
perturbation theory (considering J as a small parameter), which yields (after subtracting 
J) the eigenvalues : 



/i-|sin 2 ! -Jcos 2 f -/i-^sin 2 ! 



(10) 



By matching the results in Eq. ^|to those in Eq. [T(| in the h <^iU limit (ignoring 0(^r)) 
we can make the identification : 

At 2 

(id 

However, when the field is increased, the matching is not perfect anymore, and the effective 
exchange parameter for the ground state is reduced to 

At 2 

J -^i < 12 > 

Thus the appearence of polarons decreases the effective anti-ferromagnetic exchange between 
the carriers. 



IV. GENERALIZED HUBBARD MODEL 
A. 2-level Hubbard model with random fields 

One can improve this analysis by considering a more realistic model. The next simplest 
case is to consider two energy levels 1 and 2 (energies e\ and on each atom, and to 
allow hopping 1 «-» 1,1 <-» 2 and 2 «-> 2 between sites. Again we consider arbitrary fields h a 
and hf,. The Hamiltonian is still given by Eq. |] with the summation for % going over 1, 2. 
The number of states increases dramatically : we are dealing now with a 28 x 28 matrix 
( 28 = 6 x (3 + 1) + 4 x 1 since there are 6 pairs of different spatial states which each can 
have spin or 1, and 4 pairs of identical states which can only have spin 0). We need to 
concentrate on the lowest energy states only, treating the rest perturbatively. We ignore the 
2^2 hopping, since it affects the lowest eigenvalues only to higher order in PT. We are 
using again the simplifying assumption h a = hb = h. 

The lowest energy subspace can be identified as being spanned by {al | 61 t , al t 61 j, 
, al I bl |, al I 61 J,}. By applying second order degenerate PT in this subspace, we obtain 
the following expressions for the eigenvalues : 



S 



26! 
26! 

2ei 



2t 2 2 cos 2 § 
e2-ei+^12 
2tf 2 cos 2 § 



2^? 2 sin 2 f 
e 2 -ei+t/i2- 



2tL8in 2 f 
t/ll-ft 

2 cos^ ^ t 2 2 sin 2 ^ ^ sin 2 e 

£2— £1+J7l2 £2— ei+J7i2 — h e2—e\+Ui2+h 
l cos 2 f _ 2t 2 2 cos 2 | _ i 2 2 sin 2 | 
fla e2-ei+t^l2 £2— €1+I7l2 — /l 

^ 2t 2 2 cos 2 § 2t 2 2 sin 2 § 2t 2 j sin 2 

£2— ei+fl2 €2—ei+Ui2+h Un+h 



\ 



Wl2 

^! cos 2 f 



sin 2 f * 2 2^ 2 | 
L+C^12 — /l €2-ei+Ul2- 

2 f 2 2 sin 2 f 



"12 ""* 2 
e2-ei + !7l2+fe 

2t 2 j sin 2 S 



(13) 



/ 



One can see that in the limit ti2 — > we obtain the same results as in the 1-level Hubbard 
model analyzed in the beginning. This is a good consistency check. 

B. Magnetic properties of the ground state 

By applying second order perturbation theory to the 2-level Hubbard model (and con- 
sidering h 3> t 2 /U), we therefore obtain the following expression for the ground-state energy 



E GS = 2e 1 -h - 2tl2C °° 2 



2 e 



(14) 



2t 2 sin 2 



2t^ sin 2 i 
Uu+h' 



The angle 9 between the two fields was regarded up to this point as an external parameter. 
All the calculations so far were done under the assumption that the magnetic field was fixed 
externally. We must take however into account the fact that the field is generated by the 
actual polaron, and that although the magnitude of the field is fixed by the size of the 
polaron, the direction is free to change. Therefore, when T — > OK, O takes the value that 
minimizes the energy. Since Egs(®) = const. + A(h) sin 2 (^), with 



A(h) 



9t 2 



2t 2 



12 



2t 2 



n 



e 2 - ei + Un e 2 - ei + U 12 + h Uu + h 
the two values that minimize the value of the energy are = and = it, depending of the 
sign of the factor A(h). We can regard A, which represents the energy difference between the 
ferromagnetic and the anti-ferromagnetic configurations, as an effective exchange constant. 
For small values of the polaron field it is the anti- ferromagnetic state that dominates, whereas 
if we increase the polaron field the ground state of the system becomes ferromagnetic. 
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In order to get an idea of what parameters are essential for the transition, let us solve 
A(h) = 0, which is just a quadratic equation. The critical field is given by the only acceptable 
(positive) root: 

h c = U n (15) 

where we have defined a = (fjj) 2 and j3 = Ul2 ^~ ei ■ I n the one-level limit, thus, the 
transition disappears (we have a — > oo =>> h — > oo). The ferromagnetic configuration 
happens then as a consequence of the local magnetic fields and of the hopping to excited 
states. As a side remark we observe that while (3 depends entirely on the type of dopant used, 
a depends both on the type of dopant and on the dopant density (since all the transition 
probabilities depend on the separation of the atoms). Thus we may have a transition from 
a polaron ferromagnetic to an anti-ferromagnetic ground state configuration as the density 
is varied. It turns out that for the 2-level (hydrogenic ground state and any excited state) 
approximation, a can take virtually any value from a = for separation, to a — > oo when 
the separation becomes infinite. 



C. Many- level Hubbard Model 

One can generalize the results presented above about the ground-state of a 2-level Hub- 
bard model to the case where any number of bound excited states exists. Consider f ) 
and \ipi I) (with energies ei) to be the 1-electron states, indexed by the subscript i (i — 
for the ground-state). We consider the ground-state to be non-degenerate, however allowing 
arbitrary degeneracies for all the other states. The ground-state of the two-electron system 
in a magnetic field is then : aipo j bipo j, with energy (in the 0-th order) 2eo — h. We allow 
the hopping of an electron from the ground state of one atom to any state of the other 
atom, the coupling constants between aipo and bipi being given by toi- Expression ^ gives 
the right Hamiltonian with the summation for a being over a, b, corresponding to only two 
interacting polarons. 
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The couplings of the ground-state atpo I bipo j to various other states are given in Table 
|. Applying second order non-degenerate PT, we obtain for the ground state energy : 

E GS = 2e -h- - (Ei ~%r) cos2 f - (S* 2 4fr-rj:) ™ 2 f (16) 

We can apply the same kind of analysis as for the two-level case above, however A(h), the 
effective exchange constant whose sign dictates the magnetic configuration now becomes: 

oj-2 oy.2 oj-2 

,■ Q — Co + t>0i j e i — e + ^Oi + « ^00 + 
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TABLES 



State 


Coupling 


Energy diff. 


a{b)t/j i a(b)tp | 


±it o sin(f ) 


Uoo + h 


a(b)tpi t a(b)ip j, a(b)tpi [ a(b)ip | 


±^ i sin(f ) 


ei - e + f/oi + ^ 


a(b)A i a(b)tp j, a(b)ipi | a(b)tp | 


±%cos(|) 


- eo + f^oi 



TABLE I. The couplings of the ground state of the many-level Hubbard model to various 
excited states 



The value of the critical field above which the ferromagnetic configuration becomes ener- 
getically favorable is again given by the equation A(h c ) = 0, however this cannot be solved 
analytically anymore. 
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FIGURES 
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FIG. 1. The effective exchange A as a function of the dopant separation d. The polaron field 
was taken as h = 0.3Ry 
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FIG. 2. The critical field h c (in Ry) vs. the dopant separation d (in a#). The Log-Linear plot 

reveals the exponential behaviour of the critical field on the dopant separation. 



D. Application to Hydrogenic Centers 



We can understand the details of this change of magnetic configuration better if we derive 
the actual couplings and Coulomb terms from a simplified model of the dopant. As discussed 
in the introduction, one can use simple Hydrogenic models : two H atoms separated by a 
distance p in magnetic fields h a , h b . We obtained the 2-center overlap integrals (t 0i ) by 
using explicit machine-readable formulas that have been constructed by applying symbolic 
calculation to the (^-function method of Barnett and Coulson. The mathematical formalism 
has been described originally hereS. The symbolic calculations are described herein and 
the work cited therein. The evaluation of 1-center Coulomb integrals (Uoi) was done following 
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an approach described in some textbooks^. 

We considered Hydrogenic centers, and included transitions and Coulomb interaction 
energies between the ground state Is and the states Is, 2s, 2p z , 3s, 3p z , and 3d. Figure 
[L] plots the effective exchange A as a function of separation d between the two centers, for 
a fixed value of the polaron effective field, h = 0.3Ry, which is still small compared to the 
Rydberg. It can be seen that the effective exchange A becomes positive at a certain dopant 
separation d and thus the favorable configuration becomes ferromagnetic. Figure |] plots 
the critical field h c as a function of the dopant separation d. For d > 4cib, which is true for 
typical experimental doping densities, the minimal value of the polaron field that will provide 
a ferromagnetic interaction becomes reasonable (a few tenth of a Rydberg) and therefore we 
can conclude that our model predicts ferromagnetic interractions between polarons in DMS. 

V. SPIN HAMILTONIAN FOR MODERATE FIELDS 

Coming back to our fixed magnetic field model, we note that both the 1- and 2-level 
Hubbard models agree with the Hamiltonian of Eq. || containing the standard Heisenberg 
exchange and Zeeman terms when contributions of order 0{j^) are ignored. However when 
those contributions are taken into account, the Eq. |8| does not provide the right solutions 
anymore. The question to be asked is whether it is possible to modify the Hamiltonian so 
as to have agreement up to O(^l). It turns out that this is indeed possible. 

If we expand the terms in Eq. [7] we obtain the energies : 

( 2e + / i -^l-^^i 2e 2e-^l 2e - h - + ) (18) 

We need to add some small correction to Eq. |^ that is linear in the fields and reproduces 
the above structure. There are several ways of doing this, the simplest being to add a term 
of the form h a • Sb + hb ■ s a or h • s a x Sb or h a ■ s a + hb • Sb (or any linear combination 
thereof). It turns out that the Hamiltonian 

H H 2 = (1 - a) (ha ■ s a + h b ■ s b ) + Js a • s b + a(h a ■ s b + h b ■ s a ) (19) 
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reproduces the right structure. In the basis {a | & |, a | & I, a I & t, a I 1} it becomes : 

/ -(1 -a)h+ (i - ah) cos& -(±-f)sm& (±-f)sm& -|(l-cos0)^ 
-(i-f)smQ -^cos© + cos 0) -(f + ^)sin6 

(f-^)sin© ^(l + cos0) -^cosG (| + ^f)sin0 

y-f(l-cos0) -{{ + f)sin© (| + ^) s in0 (1 - a)h + + ah) cosG j 

After doing the necessary 45 degree rotation in the {a "J* b I, a I b f} subspace the matrix 
becomes suitable for PT and it yields the eigenvalues (to first order in J and a and after 
subtracting a common J): 

h-l sin 2 § - 2ah sin 2 § -J cos 2 f -h - { sin 2 f + 2ah sin 2 § ) (20) 

By matching the results in Eq. ^ to Eq. [18] we obtain for the parameters of the modified 
Hamiltonian : 

J=fa = £ (21) 

One can also expand the result for the 2-level Hubbard model to get a better estimate for 
the parameters. In that case one obtains : 

j=f± a = ^ + \ (22) 



VI. CONCLUSION 

The calculations presented above lead us to two conclusions : 

I They confirm once more the fact that the polarons formed in Dilute Magnetic Semi- 
conductors can interact ferromagnetically for certain dopant densities and types. This 
extends the results obtained byiS in the limit that the polarons have an important 
overlap to the situation where the two polarons do not overlap at all. Thus, the two 
qualitatively distinct effects combine in order to generate an effective ferromagnetic 
interaction of the bound magnetic polarons in a DMS. 
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II For the case where two 1-electron atoms are placed in fixed, but nonparallel external 
magnetic fields, the standard model with Heisenberg exchange and Zeeman terms is 
not a suitable approximation. Instead, the calculations above give a correction, with 
which we are able to reproduce the correct spectrum in the moderately high field 
domain (t <C h C [/). The correction represents an effective mixing of the fields at 
the two sites, which can be intuitively understood as a "transfer" of the field from 
one site to the other by the hopping electron. In the high field domain (h < U), this 
correction is not valid anymore and the correct model is a Heisenberg Hamiltonian 
with a field-dependent exchange constant. 

We will return now to Conclusion I with a few remarks. The change from antiferro- 
magnetic to ferromagnetic effective coupling in the presence of strong local fields in the 
generalized Hubbard model naturally leads one to the question whether this would actu- 
ally occur in an exact calculation. We believe it does, though the parameter values for the 
change are likely to be different. To explain our "belief", we consider the case of the hydro- 
gen molecule problem in zero field, where the issue of the effective Heisenberg Hamiltonian 
has been thoroughly discussed0'0 'i'S 
j 
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FIG. 3. The effective exchange parameter as a function of reduced distance (r/os) for 
Herring-Flicker (solid), Kolos-Wolniewicz (dots), Heitler-London (dashed) and Hubbard model 
(dot-dashed) . 

In Fig. ^| we plot the effective exchange parameter as a function of d — (r/ae) in the range 
1-8, as calculated in four different ways. The solid line represents the Herring-Flickeril (HF) 
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result, Jjjp{r ) = 1.636(r/a B ) 5 / 2 exp(— 2r /as) Ry, which is asymptotically exact (in the sense 
that Jhf{x) I ' Jexactir) —> 1 as r/ae — > oo ), while the dots are the numerically converged 
results of Kolos and WolniewiczEl (KW). Both show that J(r) is positive (antiferromagnetic) 
for all r. The dashed line is the result of the Heitler-London0 (HL) approximation, which, 
though clearly not exact, works reasonably well on this logarithmic plot for the range shown. 
It should be noted, however, that while the HL result has the right sign of J(r) for the range 
of r/a B shown, it incorrectly predicts a negative J(r) at large r/ae because it does not take 
into account polarization corrections to the ground state hydrogenic wavefunction. Finally, 
as the dot-dashed curve, is the standard Hubbard model result, J(r) = 4t 2 (r)/U, with 
t(r) = 2(1 + r/afl) exp(— r/ae) Ry and U = 5/4 Ry , calculated within the ground state 
approximation for the hydrogen wavefunctions (the generalized Hubbard model would give 
the same result in this case without external fields, as adding excited states does not alter 
the second order splitting between the lowest singlet and triplet states). 

As can be clearly seen, the Hubbard approximation overestimates J(r) by a large factor 
(this qualitative fact does not change with more refined estimates for t(r) and U ). The 
reason for the larger exchange is that Hubbard, and Hubbard-like approximations, consider 
only the kinetic exchange (due to the hopping process), which is antiferromagnetic, and 
neglect coulomb exchange which tends to favor ferromagnetism. (Such a split is often used 
in literature^) . The latter is included in the "exact" HF/KW treatment, as well as the HL 
approximation, resulting in a much lower value net (antiferromagnetic) exchange. 

We expect that inclusion of local magnetic fields hi of equal magnitude (h), which couple 
only to the electron spin is properly captured on a qualitative/semi-quantitative level by the 
extended Hubbard model. Therefore inclusion of such fields in a more accurate model will 
also result in a movement of the kinetic exchange towards ferromagnetism; consequently, 
the overall exchange will change over to ferromagnetic at some value of h. If we just add a 
field independent (ferromagnetic) coulomb term to the kinetic exchange of the generalized 
Hubbard model, the change from antiferromagnetism to ferromagnetism would be expected 
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at lower values of h than we have calculated, and make the effect we consider more relevant 
for DMS systems. 

We caution, though, that for the purely hydrogenic problem with local fields, one must 
take into account the effect of the magnetic field on the orbital wavefunction as well, and 
that will certainly affect the results, at least quantitatively. In the case of DMS, the local 
fields represent exchange fields due to interaction of hydrogenic states with local atomic 
states of the magnetic ion (Mn), and therefore their orbital effect is not the same as that 
of external magnetic fields in the H2 problem. Nevertheless, we expect these to have some 
effect on the orbital part of the wavefunction of the hydrogenic impurity0, which would 
have at least a quantitative effect on our results. 
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